Dynamics of driven vortex-antivortex matter in superconducting films with a 

magnetic dipole array 
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We investigate theoretically vortex-antivortex (v-av) matter moving in thin superconducting films 
with a regular array of in-plane magnetic dipoles. Our model considers v-av pair creation induced by 
the local current density generated by the magnetic texture and the transport current and simulates 
the dynamics of vortices and antivortices by numerical integration of the Langevin equation of 
motion. Calculations of the transport properties at zero applied field show a strong dependence 
of the v-av dynamics on the current intensity and direction. The dynamics of the v-av matter is 
characterized by a series of creation and annihilation processes, which refiect on the time dependence 
of the electrical field, and by guided motion, resulting in a zero-field transverse resistance. 
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A small magnetic dot placed in the vicinity of a thin 
superconducting film has the ability of creating vortex- 
antivortex (v-av) pairs in the superconductor and keep 
them from annihilating each otheriii^ Thanks to this, 
phases of vortex matter in which vortices and antivor- 
tices coexist can be stabilized in bilayers made of thin 
superconducting films and arrays of nanomagnets, even 
in the absence of macroscopic magnetic fields. This topic 
has attracted a great deal of attention lately. Several cal- 
culations using the Ginzbug-Landau approach have been 
carried out for arrays of polarized magnetic dots on top 
of a superconducting film, predicting a great variety of 
equilibrium phases of vortices and anti-vorticesi^ 

Most early experimental work on superconductor- 
nanomagnet hybrids were aimed at the flux pinning prop- 
erties of the magnetic dotSi^^'^ Recently, experiments on 
superconducting films with polarized magnetic dot arrays 
corroborate the existence of vortex-antivortex phases in 
these structures and show that such phases induce novel 
phenomena such as the field-induced superconductivity 
and the zero-magnetic-field ratchet efFect^i^ 

A problem that has received much less attention is 
the existence in these systems of dynamical states of 
vortex matter in which vortices and antivortices coex- 
ist. These states may result when a transport current 
is applied to the film, and the v-av matter moves. In 
this paper we study this problem theoretically. We ex- 
pect that during the motion v-av pairs annihilate and 
that new pairs are created. The fundamental question 
is how these processes, together with vortex- vortex and 
vortex-nanomagnet interactions, determine the proper- 
ties of the moving vortex matter. Here we address this 
point by using numerical simulations of a simple model 
for vortices interacting with a periodic array of perma- 
nent point dipoles. 

The details of our model is as follows. The bilayer com- 
prises a thin superconducting film of thickness d, pene- 
tration depth A d, and coherence length ^ < d. A 
square array of magnetic dipoles, with magnetic moment 



m and unit cell of dimensions ap x Up (Fig. [T]), is placed 
above the film, at a distance zq > d/2 from it. Exper- 
imentally, the space between the superconducting film 
and the dipole layer is usually filled with an insulating 
buffer layer to prevent proximity effects. Here we assume 
that the dipoles are polarized in plane, parallel to the x 
direction (m = rnx). The total energy of the vortex 
matter is given by 
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where q is the vorticity {q — I for vortices, q = —1 
for antivortices). Here we assume that no external mag- 
netic field is applied. In other words, the vortex-matter 
"neutrality" condition, qj = 0, must be satisfied. 
Uvm{i^) is the usual interaction energy between a vor- 
tex located at r and the dipole array calculated in the 
London limitj^Si The interaction energy Uwi^ij) between 
an arbitrary pair of vortices at positions and Vj is 
calculated following Clem's variational approach for the 
order parameter (r^ — ri ~ r^).— In the limit A 3> c? 
considered here, this gives J7vv(i") = eoc^ln(|ryp -I- 2^^), 
where eo = (Pq/Att^qX^ and is a variational parameter 
(in general, S.(T)). The vortex core energy can be 
calculated straightforwardly by using Clem's trial order 
parameter to compute the gain in condensation energy 




FIG. 1: (Color online) (a) Schematic view of the in-plane 
dipole array, (b) Definition of the transport current J and 
driving force F orientations. 
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due to an isolated vortex, which results in Ec = 0.375eo. 
—Ec acts as a chemical potential of the vortex system 
and thus plays an important role in the nucleation of 
vortex-antivortex pairs. 

Our numerical algorithm is based on the assumption 
that the time evolution of the vortex matter can be bro- 
ken into two distinct parts: motion of the vortices and 
antivortices, and creation of v-av pairs. The first is gov- 
erned by the Langevin equations 
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dt 
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where rj is the Bardeen-Stephen friction coefficient, j 
runs over all vortices and anti-vortices present at time 
t, Fj = {qj(pQ/c)3 X z is the driving force, and Tj is the 
random force appropriate to temperature T. J is the 
applied current density which makes and angle a with 
the positive y-axis, as shown in Fig. [TJb). These equa- 
tions are solved numerically by a finite difference method. 
Eventually, a vortex and an antivortex collide and anni- 
hilate each other. This is taken into account by assuming 
that such a collision occurs every time the vortex and the 
antivortex are separated by a distance ^c- 

The procedure for the second part, v-av pair creation, 
is to freeze the vortex matter configuration after every 
Langevin dynamics time step and create v-av pairs using 
a Monte Carlo procedure. This assumes that the actual 
nucleation process takes place in a time interval so short 
that the motion of the vortex matter during it can be 
neglected. Trial v-av pairs are then placed at every point 
of an auxiliary grid, with unit cell of dimensions ac x ac 
(oc — £,c), commensurate with the dipole array. The vor- 
tex and antivortex are placed at a distance from each 
other at a given axis, with the pair center of mass fixed at 
the grid point. At high T, this axis is chosen randomly, 
whereas at low T it is determined by the direction of 
the total current density calculated at the grid point (at 
that instant of time) as to favor pair unbinding. Then, a 
Metropolis code is used to either accept or reject the pair. 
The energy entering this algorithm includes the energy 
due to the interaction with the transport current, that is 
Et = E — J2j Fj • . Such a procedure accounts for both 
deterministic (induced by the local current distribution) 
and thermal pair nucleationii^ 

The simulations are carried out on a L x L section of 
the bilayer, with L = 8ap, and assuming periodic bound- 
ary conditions in both x and y directions. The time 
step used in the numerical integration is dt — 2 ■ lO^^^o 
{to = (jiQeQ / (rjap)) and we chose zq = 0.2ap. Equilib- 
rium configurations of the vortex matter for small T and 
J = 0, were obtained minimizing E (Eq. [1]) by means 
of a simulated annealing scheme. Fig. [5] shows the equi- 
librium states of the superconducting film obtained by 
this method. The bottom panel presents the number of 
vortex-antivortex pairs per dipole (ripair) as a function of 
the magnetic moment of the dipoles for several values of 
^c- These curves are characterized by plateaus at integer 
Upair- Examples of stable vortex-antivortex configura- 
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FIG. 2; (Color Online) (Top) Vortex(black)-antivortex(red) 
configurations for (a) m = 1.0(^0^0 (n-pair ~ 1), (b) m = 
1.65<^o2o {ripair = 2) and (c) m = 2.250o2o (npair = 3). (Bot- 
tom) Number of v-av pairs per dipole, ripair, versus magnetic 
moment of the dipoles, m, for different values of Ac- 



tions for £,c — 7x 10~^ap are depicted in the upper panels. 
Notice that by increasing (for instance, by increasing 
the temperature) the magnetic moment necessary for the 
creation of a new v-av pair per dipole decreases. This is 
consistent with the fact that the force necessary to sep- 
arate a v-av pair decreases with ^c- Therefore, for lower 

values, a strong magnetic moment is necessary to hold 
the pair from annihilating each other. A similar temper- 
ature dependence of the critical magnetic moment has 
been also observed in numerical calculation for magnetic 
dot arrays in the off-plane geometry^. 

Now we discuss the response of the vortex matter to 
an applied current flow. For simplicity, we shall consider 
hereafter only the cases where each dipole stabilize one 
v-av pair in equilibrium, by assuming — 7 x lO^^Up 
and m — 1.0(/)ozo- Fig-ISKa) presents the x and y compo- 
nents of the mean electric field induced by vortex motion, 
expressed by the time average of E = 0o 

Q^Vi x as a 

function of the current density J (in units of the depairing 
current Jd) for different orientations of J. Averages were 
taken over 15 • 10^ time steps after a relaxation waiting 
time of 5 • 10^ time steps. The Ex{J) and Ey{J) curves 
evidence a strong dependence of the v-av dynamics on the 
current orientation a. Such anisotropy becomes clearer 
when one plots the critical current Jc as a function of a. 
Vortices are depinned more easily when J is at an angle 
a — 180°. At this direction, 3 ~ Jy induces a Lorentz 
force on the vortices (antivortices) parallel to +x {—x) 
in such a way as to separate the v-av pairs. Such separa- 
tion is favored by the supercurrent generated right below 
each dipole, which is also parallel to —y and is respon- 
sible for the stabilization of the v-av pairs. At a = 90°, 
on the other hand, the applied current has to overcome 
alone the v-av mutual attraction, with no help from the 
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FIG. 3: (Color online) Time averages of the x and y compo- 
nents of the electric field (top) and of the v-av pair density 
(Bottom) as functions of the current density applied at differ- 
ent orientations a. Inset: angular dependence of the critical 
current, Jc{ce). 



dipole- induced supercurrent. For a < 90°, the compo- 
nent Jy inverts sign and tends to bind the v-av pair. At 
a — 75°, Jx is still strong enough to unbind the pairs. 
At lower a, however, all v-av pairs annihilate each other 
and the current necessary to create new pairs and stab- 
lish a dynamical state of moving v-av matter is close to 
Jd- At such current range, our model is no longer valid 
and therefore Jc could not be estimated for a < 75°. 
Fig. El (b) presents the mean v-av pair density {upair) as 
a function of J. Interestingly, for all angles studied the 
onset of v-av matter motion is accompanied by an in- 
crease in (upair) with respect to the static (equilibrium) 
value (npair) = 1- 

Here we analyze in more detail the dynamics of vor- 
tices and antivortices in the film and how it reflects on 
the macroscopic electrical response and the mean pair 
density. We shall restrict ourselves to four directions, 
a = 180°, a = 150°, a = 120° and a = 90°. Fig. g] 
present the time evolution of the number of v-av pairs 
per dipole and the x and y components of the mean elec- 
tric field for fixed J values. For J > Jc(a), we observed 
steady states of the moving v-av matter characterized 
by an oscillatory time dependence of and Ey, even 
though the applied current is constant. For a — 90°, 
npair(t) rapidly reaches 1.5 and keeps essentially con- 
stant at this value. The corresponding steady state is 
illustrated in Fig. [5]^ a), where we present contour plots 
of the time-averaged vortex-antivortex density, defined 
as p(r) = {J2i 9i<5(r — ri{t)))t- The average is taken at a 
fixed J over iVg time steps. The trajectories of vortices 
(antivortices) are given by the light (dark) color shades. 
As it is clear, vortices and antivortices follow distinct 
tracks, perpendicular to the applied current, and collide 
very rarely. 
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FIG. 4: Time series of the electric field Ex, Ey and v-av pair 
density Upa,,- for a = 90° (a), 120° (b), 150° (c) and 180° (d). 



For a > 90°, the dynamical steady states are much 
reacher. npair{t) is also oscillatory, indicating a cyclic 
series of creation and annihilation of v-av pairs. The 
general picture is as follows: when a current J > Jc{a) 
is applied, the vortices and antivortices stabilized by the 
dipoles in the equilibrium (J = 0) state, which we shall 
refer to as first-generation vortices, depin and move apart 
from each other, leaving space for the creation of new 
(second-generation) v-av pairs right below the dipoles. 
At this point, npair{t) reaches twice its original value 
{npair{t) = 2). These new pairs, give an additional push 
to the preexistent ones, inducing a sudden increase in 
E|. All first-generation vortices and antivortices run to- 
wards annihilation, which will take place at the midpoint 
between first-neighbor dipoles, for 135° < a < 180°, 
or at the face center of the dipole array unit cells, for 
105° < a < 135°, as revealed by the p(r) plots pre- 
sented in Fig. [Sl^a). Just before annihilation, however, 
new v-av pairs spawn at the above mentioned positions 
and subsequently annihilate the first-generation vortices 
and antivortices. We believe these extra pairs arise from 
the fact that our model assumes a rigid vortex core. At 
distances of order ^, the v-av annihilation in supercon- 
ductors is characterized by strong vortex core deforma- 
tion. In our model, the extra v-av pairs appear naturally 
as a way of compensating the vortex core rigidity and 
accelerate the annihilation process. As a side effect, a 
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FIG. 5: (Color online) (a) Time- averaged vortex density for 
J — 0.5Jd and different current orientations a. Light (dark) 
shades correspond to vortices (antivortices). The circles in- 
dicate creation (★) and annihilation (f) of v-av pairs, (b) 
Electric field direction as a function of a for several driv- 
ing current intensities (symbols). The — a line is shown 
for comparison. The cartoons illustrate the dominating v-av 
motion for 9 ~ 180° and 61 ~ 90°. 

spurious, high-frequency noise arises in Ex{t) and Ey{t). 
Notwithstanding, we have checked that such noise has 
only a neglectable influence on the macroscopic quantities 
shown in Fig. [S] Finally, after the annihilation process, 
the second-generation v-av pairs become first-generation 
ones and the whole process repeats. Notice that the life- 



time, ^L, for a vortex of the first or second generation is 
always larger than the oscillation period P (see Fig. 2]). 

The angular dependence of the v-av dynamics de- 
scribed above suggests that motion occurs preferentially 
along the x directions for a close to 180°. Such guided 
motion of the v-av matter should reflect directly on the 
electric field E generated by vortex motion, giving rise 
to a deflection of E with respect to J and thus inducing 
transverse resistance. To demonstrate this statement, we 
plot the direction of E, 9, as a function of the direction 
of J, a, for several values of J (left panel of Fig. [5]). In- 
deed, 9{a) diverges considerably from the = a line. 
For a close to 180°, 9 tends to lock at 6* ~ 180° in the 
whole J range, whereas for a = 75°, E locks at = 90°. 
This transverse resistance resulting from the guided mo- 
tion of v-av matter is essentially different from that ob- 
served in superconducting films with non-magnetic pin- 
ning arraysii. and other Hall-like effects in the sense that 
here no macroscopic magnetic field is required. 

In summary we have studied the dynamics of vor- 
tices and antivortices in superconducting films interact- 
ing with an array of in-plane magnetic dipoles by means 
of a model which explicitly takes into account creation 
and annihilation of v-av pairs. The current-driven dy- 
namics of the vortex-antivortex matter at zero applied 
field is characterized by a series of creation and annihi- 
lation processes producing an oscillatory behavior of the 
time dependent electric field. Other predictions of our 
model include the anisotropy in the critical current and 
the transverse electric field induced by vortex-antivortex 
guidance. These results can be checked out experimen- 
tally by means of conventional transport measurements 
at zero applied magnetic field or by direct high-resolution 
magneto-optical imaging. We anticipate that the zero- 
field guidance effect should also occur for other dipole 
polarizations, such as the off-plane geometry, as long as 
the magnetic texture captures v-av pairs. 
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